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ABSTRACT 

Context. Detecting massive satellites that orbit extrasolar planets has now become feasible, which led naturally to questions about 
the habitability of exomoons. In a previous study we presented constraints on the habitability of moons from stellar and planetary 
| illumination as well as from tidal heating. 

Aims. Here I refine our model by including the effect of eclipses on the orbit-averaged illumination. I then apply an analytic approxi- 
mation for the Hill stability of a satellite to identify the range of stellar and planetary masses in which moons can be habitable. Moons 
in low-mass stellar systems must orbit their planet very closely to remain bounded, which puts them at risk of strong tidal heating. 
Methods. I first describe the effect of eclipses on the stellar illumination of satellites. Then I calculate the orbit-averaged energy 
flux, which includes illumination from the planet and tidal heating to parametrize exomoon habitability as a function of stellar mass, 
planetary mass, and planet-moon orbital eccentricity. The habitability limit is defined by a scaling relation at which a moon loses its 
water by the runaway greenhouse process. As a working hypothesis, orbital stability is assumed if the moon's orbital period is less 
' than 1/9 of the planet's orbital period. 

, Results. Due to eclipses, a satellite in a close orbit can experience a reduction in orbit-averaged stellar flux by up to about 6 %. The 

smaller the semi-major axis and the lower the inclination of the moon's orbit, the stronger the reduction. I find a lower mass limit of 
i-C 1 a 0.2 M Q for exomoon host stars that allows a moon to receive an orbit-averaged stellar flux comparable to the Earth's, with which 

it can also avoid the runaway greenhouse effect. Precise estimates depend on the satellite's orbital eccentricity. Deleterious effects on 
exomoon habitability may occur up to « 0.5 M e if the satellite's eccentricity is > 0.05. 

Conclusions. Although the traditional habitable zone lies close to low-mass stars, which allows for many transits of planet-moon 
£jf ■ binaries within a given observation cycle, resources should not be spent to trace habitable satellites around them. Gravitational pertur- 

bations by the close star, another planet, or another satellite induce eccentricities that likely make any moon uninhabitable. Estimates 
for individual systems require dynamical simulations that include perturbations among all bodies and tidal heating in the satellite. 
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O; 

O ! The detection of dozens of Super-Earths and Jupiter-mass plan- 2 - 1 ■ Orbit-averaged energy flux 

Q\ ; ets in the stellar habitable zone naturally makes us wonder Xo assess the habitability of a satellite, I estimated the global 

about the habitability of their moons (Williams et al. 1997; average 
Kaltenegger 2010). So far, no extrasolar moon has been con- 

~~ < firmed, but dedicated surveys are underway (Kipping et al. ^ (I -a) I 7rR 2 a p 

. 2012). Several studies have addressed orbital stability of extra- Ff° = s ' ■ 

'km ' so ^ ar satellite systems (Donnison 2010; Weidner & Horne 2010) 167raJ p Jl - ej, 
y\ and tidal heating has been shown to be an important energy 

H source in satellites (Reynolds et al. 1987; Scharf 2006; Cassidy ^ p o"SB(r p q ) (l - Ss ) 

. ■ et al. 2009). In a recent study (Heller & Barnes 2012, HB12 + 4 + " s (1) 

in the following), we have extended these concepts to the illu- ps 

mination from the planet, i.e. stellar reflected light and thermal of its energy flux over one stellar orbit; where u is stellar lumi . 
emission, and presented a model that invokes the runaway green- nosityi a ^ the sem i- ma jor axis of the planet's orbit about the star, 
house effect to constrain exomoon habitability. flps the sem i_ ma j or axis of the satellite's orbit about the planet, 
Earth-mass moons about Jupiter-mass planets have been e„ P the circumstellar orbital eccentricity, R p the planetary radius, 
shown to be dynamically stable for the lifetime of the solar sys- a p and ff s are the albedos of the planet and the satellite, respec- 
tem in systems where the stellar mass is greater than > 0.15 M tively, 7^ is the planet's thermal equilibrium temperature, h s 
(Barnes & O'Brien 2002). I consider here whether such moons the satellite's surface-averaged tidal heating flux, cr SB the Stefan- 
could actually be habitable. Therefore, I provide the first study Boltzmann constant, and x s is the fraction of the satellite's orbit 
of exomoon habitability that combines effects of stellar illumi- that is not spent in the shadow of the planet. Tidal heating /j s 
nation, reflected stellar light from the planet, the planet's thermal depends on the satellite's eccentricity e ps , on its radius R s , and 
emission, eclipses, and tidal heating with constraints by orbital strongly on a ps . I avoid repeating the various approaches avail- 
stability and from the runaway greenhouse effect. able for a parametrization of /i s and refer the reader to HB12 
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and Heller et al. (2011), where we discussed the constant-time- 
lag model by Leconte et al. (2010) and the constant-phase-lag 
model from Ferraz-Mello et al. (2008). For this study, I arbitrar- 
ily chose the constant-time-lag model with a tidal time lag of 
the model satellites similar to that of the Earth, i.e. r s = 638 s 
(Neron de Surgy & Laskar 1997). 

In HB12 we derived Eq. (1) for x s = 1, thereby neglect- 
ing the effect of eclipses on the average flux. Here, I explore 
the decrease of the average stellar flux on the satellite due to 
eclipses. Eclipses occur most frequently - and will thus have 
the strongest effect on the moon's climate - when the satellite's 
circum-planetary orbit is coplanar with the circumstellar orbit. If 
the two orbits are also circular and eclipses are total 



x s = 1 - R p /(na ps ) . 



(2) 



Applying the Roche criterion for a fluid-like body (Weidner & 
Home 2010), I derive 
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2.44tt R, \ M, 



< x s < 1 



(3) 



where x s = 1 can occur when a moon's line of nodes never 
crosses the planet's disk, i.e. in wide orbits. For an Earth-sized 
satellite about a Jupiter-sized planet Eq. (3) yields x s > 79 %. 
Due to tidal effects, the moon's semi-major axis will typically 
be > 5R p , then x s > 1 - l/(57r) « 93.6 %. In other words, 
eclipses will cause an orbit-averaged decrease in direct stellar 
illumination of » 6.4 % at most in realistic scenarios. Using 
a constraint similar to Eq. (3), Scharf (2006) derived eclipsing 
effects of the same order of magnitude. 

Eclipses will be most relevant for exomoons in low-mass 
stellar systems for two reasons. Firstly, tides raised by the star 
on the planet will cause the planet's obliquity to be eroded in 
<K 1 Gyr (Heller et al. 2011), and since moons will orbit their 
planets in the equatorial plane (Porter & Grundy 201 1), eclipses 
will then be most likely to occur. Secondly, exomoons in low- 
mass star systems must be close to their planet to ensure stabil- 
ity (see Sect. 2.3), hence, eclipses will always be total. Below, I 
compute the effect of an eclipse-induced decrease of stellar irra- 
diation for some examples. 

2.2. The circumstellar habitable zone 

The irradiation habitable zone (IHZ) is defined as the circum- 
stellar distance range in which liquid surface water can persist 
on a terrestrial planet (Dole 1964). On the one hand, if a planet 
is too close to the star, its atmosphere will become saturated with 
H2O. Photodissociation then drives an escape of hydrogen into 
space, which desiccates the planet by turning it into a runaway 
greenhouse (Kasting et al. 1993). On the other hand, if the planet 
is too far away from the star, condensation of the greenhouse gas 
CO2 will let any liquid surface water freeze (Kasting et al. 1993). 
This picture can be applied to exomoons as well. However, in- 
stead of direct stellar illumination only, the star's reflected light 
from the planet, the planet's thermal emission, and tidal heating 
in the moon need to be taken into account for the global energy 
flux (HB12). 

Analytic expressions exist that parametrize the width of the 
IHZ as a function of stellar luminosity. I used the set of equa- 
tions from Selsis et al. (2007) to compute the Sun-Earth equiva- 
lent distance Z oe between the host star and the planet - or in this 
case: the moon. To a certain extent, the constraints on exomoon 
habitability will therefore still allow for moons with massive car- 
bon dioxide atmospheres and strong greenhouse effect to exist 



in the outermost regions of the IHZ but I will not consider these 
special cases. 

Exomoon habitability can be constrained by an upper limit 
for the orbit-averaged global flux. Above a certain value, typ- 
ically around 300 W/m 2 for Earth-like bodies, a satellite will 
be subject to a runaway greenhouse effect. This value does not 
depend on the atmospheric composition, except that it contains 
water (Kasting 1988). As in HB12, 1 used the semi-analytic ex- 
pression from Pierrehumbert (2010) to compute this limit, Frg- 
For computations of the satellite's radius, I used the model of 
Fortney et al. (2007) with a 68 % rock-to-mass fraction, similar 
to the Earth. 



2.3. Orbital stability 

A body's Hill radius R^ is the distance range out to which the 
body's gravity dominates the effect on a test particle. In realis- 
tic scenarios the critical semi-major axis for a satellite to remain 
bound to its host planet is merely a fraction of the Hill radius, 
i.e. a ps < f Rh (Holman & Wiegert 1999), with a conservative 
choice for prograde satellites being / = 1/3 (Barnes & O'Brien 
2002). If the combined mass of a planet-moon binary measured 
by radial-velocity is much smaller than the stellar mass, but still 
much largerthan the satellite's mass, then P ps /P,p < 1/9, where 
P t p is the circumstellar period of the planet-moon pair and P ps 
is the satellite's orbital period about the planet (Kipping 2009). 
With M t » Mp » M s and the planet-moon barycenter orbit- 
ing at an orbital distance Z 0ffi to the star, Kepler's third law gives 



< h 



1 (M p + M s ) 



81 (M t +Mp) 



1/3 



(4) 



With this dynamical restriction in mind, we can think of sce- 
narios in which the mass of the star is very small, so that / OQ lies 
close to it, and a moon about a planet in the IHZ must be very 
close to the planet to remain gravitationally bound. Then there 
will exist a limit (in terms of minimum stellar mass) at which 
the satellite needs to be so close to the planet that its tidal heat- 
ing becomes so strong that it will initiate a greenhouse effect and 
will not be habitable. Below, I determine this minimum mass for 
stars to host potentially habitable exomoons. 

TheP ps /P„p < 1/9 constraint is used as a working hypoth- 
esis, backed up by numerical studies (Barnes & O'Brien 2002; 
Domingos et al. 2006; Cassidy et al. 2009) and observations of 
solar system satellites. But more extreme scenarios such as retro- 
grade, Super-Ganymede exomoons in a wide orbit about planets 
in the IHZ of low-mass stars are theoretically possible (Namouni 
2010). I did not consider the evolution of the planet's radius, 
its co-rotation radius, or the satellite's orbit about the planet. 
Indeed, those aspects would further reduce the likelihood of hab- 
itable satellites for given star-planet masses, because they would 
remove moons from the planet's orbit (Barnes & O'Brien 2002; 
Domingos et al. 2006; Sasaki et al. 2012). 

3. Results 

3. 1 . Eclipse-induced decrease of stellar illumination 

In Fig. 1 I show two examples for x s < 1. The ordinate gives 
the perpendicular displacement r± of the satellite from the planet 
in units of planetary radii and as a function of circumstellar or- 
bital phase tp tp . Thus, if r± < 1 and if the moon is behind the 



1 Computations were performed with my python code 
exomoon. py. Download via www.aip.de/People/RHeller. 
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Fig. 1. Perpendicular displacement of a moon as seen from the star for two scenarios. The thin and highly oscillating curve corre- 
sponds to the orbit of a satellite in a Europa-wide orbit (in units of R p ) about a Jupiter-sized planet at 1 AU from a Sun-like star. The 
satellite's orbit about the Jovian planet is inclined by i = 45°. The thicker gray line shows the Miranda-wide orbit of a satellite 
orbiting a Neptune-mass object at as 0.16 AU from a 0.4 M G star, i.e. in the center of the IHZ. Here, the moon orbits in the same 
plane as the planet orbits the low-mass host star, i.e. i = 0°. Time spent in eclipse is emphasized with a thick, black line. Both 
orbits are normalized to the circumstellar orbital phase <p tp . 



planet (and not in front of it), an eclipse occurs. Effects of the 
planet's penumbra are neglected. 

In one scenario (thin, highly oscillating line) I placed a satel- 
lite in a Europa-wide orbit, i.e. at 9.7 R p , about a Jupiter-sized 
planet in 1 AU distance from a Sun-like star. I applied a stel- 
lar orbital eccentricity of 0.3 and the moon's orbit is tilted by 
45° against the circumstellar orbit. In this specific scenario, 
the eccentricity causes the summer (0 < <p tp < 0.18 
and 0.82 < < 1) to be shorter than the winter 

(0.18 < (p* p < 0.82), and the inclination causes eclipses to 
occur only during a small part of the 365 d-period circumstel- 
lar orbit, namely around tp* p w 0.18 and <p * p a 0.82. Owing 
to the relatively wide satellite orbit, occultations are short com- 
pared to the moon's orbital period. Thus, effects of eclipses on 
the satellite's climate are small, in this case I find x s « 99.8 %. 

In the second scenario (gray line), I put a satellite in a close 
orbit at 5.1 R p , comparable to Miranda's orbit about Uranus, 
about a Neptune-mass planet, which orbits a 0.4 M star in the 
center of the IHZ at « 0. 16 AU. The planet's circumstellar orbit 
with a period of roughly 36 days is circular and the moon's three- 
day orbit is coplanar with the stellar orbit. In this case, transits 
(thick black line) occur periodically once per satellite orbit and 
effects on its climate will be significant with x s a 93.8%. 

3.2. Minimum stellar masses for habitable exomoons 

In Fig. 2 contours present limiting stellar masses (abscissae) 
and host planetary masses (ordinate) for satellites to be habit- 
able. I considered two moons: one with a mass ten times that 
of Ganymede (10 Mc a , left panel) and one with the mass of the 
Earth (M e , right panel), respectively. Variation of the satellite's 
mass changes the critical orbit-averaged flux Frq (see title of 



each panel). Contours of Ff = Frg are drawn for five dif- 
ferent eccentricities of the satellite, four of which are taken from 
the solar system: e lo = 0.0041, e Ca = 0.0074, e Eu = 0.0094, 
and eji = 0.0288 for Io, Callisto, Europa, and Titan, respec- 
tively. I also plot a threshold for a hypothetical eccentricity of 
0.05. Moons in star-planet systems with masses M„ and M p left 
to a contour for a given satellite eccentricity are uninhabitable. 
As a confidence estimate of the tidal model, dashed contours 
provide the limiting mass combinations for strongly reduced 
(t s = 63.8 s) and strongly enhanced (t s = 6380 s) tidal dissi- 
pation in the moon for the e ps = 0.05 example. 

Satellites with the lowest eccentricities can be habitable even 
in the lowest-mass stellar systems, i.e. down to 0.1 M and be- 
low for extremely circular planet-moon orbits. With increasing 
eccentricity, however, the satellite needs to be farther away from 
the planet to avoid a greenhouse effect, consequently it needs to 
be in a wider orbit about its planet, which in turn means that the 
planet needs to be farther away from the star for the moon to 
satisfy Eq. (4). Thus, the star needs to be more massive to have a 
higher luminosity and to reassure that the planet-moon binary is 
within the IHZ. This trend explains the increasing minimum stel- 
lar mass for increasing eccentricities in Fig. 2. For e ps = 0.05 
this minimum stellar mass is between 0.3 M for a 10 Mean-mass 
moon (left panel) and 0.36 M for an Earth-mass satellite (right 
panel). Uncertainties in the parametrization of tidal dissipation 
increase this limit to almost 0.5 M . 



4. Conclusions 

Eclipses of moons that are in tight orbits about their planet can 
substantially decrease the satellite's orbit-averaged stellar illu- 
mination. Equation (2) can be used to compute the reduction by 
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Fig. 2. Contours of maximum stellar masses (abscissa) for two possibly habitable moons with given planetary host masses (ordinate). 
Solid lines correspond to the satellite's time lag t s of 638 s and one out of five eccentricities, namely those of Io, Callisto, Europa, 
Titan, and 0.05, from left to right. Dashed lines refer to e ps = 0.05 and variation of t s by a factor of ten. A satellite with host star 
and host planet masses located left to its respective eccentricity contour is uninhabitable. 



total eclipses for circular and coplanar orbits. In orbits similar to 
the closest found in the solar system, this formula yields a reduc- 
tion of about 6.4 %. In orbits wider than 32 planetary radii, the 
reduction is < 1 %. In tight orbits, illumination from the planet 
partly compensates for this reduction (HB12). 

M dwarfs with masses < 0.2 M Q can hardly host habit- 
able exomoons. These moons' orbits about their planet would 
need to be almost perfectly circularized, which means that they 
would need to be the only massive moon in that satellite system. 
Moreover, the nearby star will excite substantial eccentricities 
in the moon's orbit. This increases the minimum mass of host 
stars for habitable moons. To further constrain exomoon habit- 
ability, it will be necessary to simulate the eccentricity evolution 
of satellites with a model that considers both A^-body interaction 
and tidal evolution. Such simulations are beyond the scope of 
this communication and will be conducted in a later study. 

The full Kozai cycle and tidal friction model of Porter & 
Grundy (201 1) demonstrates very fast orbital evolution of satel- 
lites that orbit M stars in their IHZ. The star forces the satel- 
lite orbit to be eccentric, thus the moon's semi-major axis will 
shrink due to the tides on timescales much shorter than one Myr. 
Eventually, its orbit will be circular and very tight. Indeed, a 
moon's circum-planetary period in the IHZ about an M0 star 
will be even shorter than the P* p /9 I used as a working hypothe- 
sis (see Fig. 3 in Porter & Grundy 201 1). 

For low-mass stars whose IHZs are close in and whose plan- 
ets' obliquities are eroded in <k 1 Gyr, eclipses will reduce the 
orbit-averaged stellar illumination of satellites by about a few 
percent. To compensate for this cooling, the planet-moon duet 
would need to be closer to the star, but then the moon would be- 
come uninhabitable by initiating a greenhouse effect. This con- 
tradiction makes my estimates for the minimum stellar mass for 
a habitable moon star more robust. 

Additionally, moons about low-mass stars will most likely 
not experience seasons because they would orbit their planet 
in its equatorial plane (Porter & Grundy 2011). The paucity 
of Jovian planets in the IHZ of low-mass stars (Borucki et al. 
201 1) further decreases the chances for habitable Earth-sized ex- 
omoons in low-mass stellar systems. 

I conclude that stars with masses < 0.2 M Q cannot host hab- 
itable moons, and stars with masses up to 0.5 M Q can be affected 



by the energy flux and orbital stability criteria combined in this 
paper. A model that couples gravitational scattering with tidal 
evolution is required to further constrain exomoon habitability 
about low-mass stars. 
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